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Abstract 

The Turaev-Viro invariants are a powerful family of topological invari¬ 
ants for distinguishing between different 3-manifolds. They are invaluable 
for mathematical software, but current algorithms to compute them re¬ 
quire exponential time. 

The invariants are parameterised by an integer r > 3. We resolve the 
question of complexity for r = 3 and r = 4, giving simple proofs that 
computing Turaev-Viro invariants for r = 3 is polynomial time, but for 
r = 4 is #P-hard. Moreover, we give an explicit fixed-parameter tractable 
algorithm for arbitrary r, and show through concrete implementation and 
experimentation that this algorithm is practical—and indeed preferable— 
to the prior state of the art for real computation. 

Keywords: Computational topology, 3-manifolds, invariants, #P-hardness, 
parameterised complexity 


1 Introduction 

In geometric topology, testing homeomorphism (topological equivalence) is a 
fundamental algorithmic problem. However, beyond dimension two it is re¬ 
markably difficult. In dimension three—the focus of this paper—an algorithm 
follows from Perelman’s proof of the geometrisation conjecture [T^], but it is ex¬ 
tremely intricate, its complexity is unknown and it has never been implemented. 

As a result, practitioners in computational topology rely on simpler topolog¬ 
ical invariants —computable properties of a topological space that can be used 
to tell different spaces apart. One of the best known invariants is homology, but 
for 3-manifolds (the 3-dimensional generalisation of surfaces) this is weak: there 
are many topologically different 3-manifolds that homology cannot distinguish. 
Therefore major software packages in 3-manifold topology rely on invariants 
that are stronger but more difficult to compute. 

In the discrete setting, among the most useful invariants for 3-manifolds are 
the Turaev- Viro invariants [S]. These are analogous to the Jones polynomial 
for knots: they derive from quantum field theory, but offer a much simpler 
combinatorial interpretation that lends itself well to algorithms and exact com¬ 
putation. They are implemented in the major software packages Regina [5] and 
the Manifold Recogniser mun], and they play a key role in developing census 
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databases, which are analogous to the well-known dictionaries of knots [Sill]. 
Their main difficulty is that they are slow to compute: current implementations 
0115] are based on backtracking searches, and require exponential time. 

The purpose of this paper is threefold: (i) to introduce the Turaev-Viro in¬ 
variants to the wider computational topology community; (ii) to understand the 
complexity of computing these invariants; and (iii) to develop new algorithms 
that are suitable for practical software. 

The Turaev-Viro invariants are parameterised by two integers r and q, with 
r > 3; we denote these invariants by "TYr,q- A typical algorithm for computing 
TVj.^g will take as input a triangulated 3-manifold, composed of n tetrahedra 
attached along their triangular faces; we use n to indicate the input size. For 
all known algorithms, the difficulty of computing TVr-^g grows significantly as r 
increases (but in contrast, the difficulty is essentially independent of q). 

Our main results are as follows. 

• Kauffman and Lins [3] state that for r = 3,4 one can compute via 

“simple and efficient methods of linear algebra”, but they give no details 
on either the algorithms or the complexity. We show here that in fact the 
situations for r = 3 and r = 4 are markedly different: computing TVr,q 
for orientable manifolds and r = 3 is polynomial time, but for r = 4 is 
#P-hard. 

• We give an explicit algorithm for computing TVi-^q for general r that 
is fixed-parameter tractable (FPT). Specifically, for any fixed r and any 
class of input triangulations whose dual graphs have bounded treewidth, 
the algorithm has running time linear in n. Furthermore, we show through 
comprehensive experimentation that this algorithm is practical —we im¬ 
plement it in the open-source software package Regina 0, run it through 
exhaustive census databases, and find that this new FPT algorithm is 
comparable to—and often significantly faster than—the prior backtrack¬ 
ing algorithm. 

• We give a new geometric interpretation of the formula for TV^^g, based on 
systems of “normal arcs” in triangles. This generalises earlier observations 
of Kauffman and Lins for r = 3 based on embedded surfaces 0 , and offers 
an interesting potential for future algorithms based on Hilbert bases. 

The ^P-hardness result for r = 4 is the first classical hardness result for 
the Turaev-Viro invariants^ However, the proofs for this and the polynomial¬ 
time r = 3 result are simple: the algorithm for r = 3 derives from a known 
homological formulation and the result for r = 4 adapts Kirby and Melvin’s 
NP-hardness proof for the more complex Witten-Reshetikhin-Turaev invariants 

m- 

The FPT algorithm for general r is significant in that it is not just the¬ 
oretical, but also practical—and indeed preferable —for real software. It was 
previously known that computing TVr,q is FPT 0, but that prior result was 
purely existential, and would lead to infeasibly large constants in the running 
time if translated to a concrete algorithm. More generally, FPT algorithms do 
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not always translate well into practical software tools, and this paper is signif¬ 
icant in giving the first demonstrably practical FPT algorithm in 3-manifold 
topology. 

2 Preliminaries 

Let M be a closed 3-manifold. A generalised triangulation of M is a collection 
of n abstract tetrahedra Ai,...,A„ equipped with affine maps that identify 
(or “glue together”) their An triangular faces in pairs, so that the underlying 
topological space is homeomorphic to M. 

In particular, as a consequence of the face identifications, it is possible that 
several vertices of the same tetrahedron may be identified together (and likewise 
for edges and triangles). Indeed, it is common in practical applications to have 
a one-vertex triangulation, in which all vertices of all tetrahedra are identified 
to a common point. In general, the 4n tetrahedron vertices are partitioned into 
equivalence classes according to how they are identified together; we refer to 
each such equivalence class as a single vertex of the triangulation, and likewise 
for edges and triangles. 

Generalised triangulations are widely used across major 3-manifold software 
packages. They are (as the name suggests) more general than simplicial com¬ 
plexes, which allows them to express a rich variety of different 3-manifolds using 
very few tetrahedra. For instance, with just n < 11 tetrahedra one can create 
13400 distinct prime orientable 3-manifolds umi. 

2.1 The Turaev-Viro invariants 

Let T be a generalised triangulation of a closed 3-manifold M, and let r and q be 
integers with r > 3, 0 < g < 2r, and gcd(r, q) = 1. We define the Turaev-Viro 
invariant TVr,g(T) as follows. 

Let V, E, F and T denote the set of vertices, edges, triangles and tetrahedra 
respectively of the triangulation A. Let / = {0,1/2,1,3/2,..., (r — 2)/2}; note 
that |/| = r — 1. We define a colouring of A to be a map 0: E —>■ I; that is, 9 
“colours” each edge of A with an element of /. A colouring 9 is admissible if, for 
each triangle of A, the three edges ei, 62 , and 63 bounding the triangle satisfy: 

• the parity condition 9{ei) 6 ^( 62 ) -I- ^(ea) G Z; 

• the triangle inequalities 9{ei) < 6 ^( 62 ) -I- ^( 63 ), 9{e2) < ^(ei) -t- 9 ( 03 ), and 
^*( 63 ) < 9(ei) -h 9 ( 62 ); and 

• the upper bound constraint 9(ei) ^( 62 ) -|- 9 ( 63 ) < r — 2. 

More generally, we refer to any triple {i,j, k) G I x I x I satisfying these three 
conditions as an admissible triple of colours. 

For each admissible colouring 9 and for each vertex v G V, edge e G E, 
triangle f G F or tetrahedron t G T, we define weights |u|e, \e\e, \f\e, |t|e € C. 

Our notation differs slightly from Turaev and Viro m-: most notably, Tu¬ 
raev and Viro do not consider triangle weights \f\g, but instead incorporate an 
additional factor of 1 / 1 ^^ into each tetrahedron weight \t\g and \t'\g for the two 
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Figure 1: Edge colours of a tetrahedron. 


tetrahedra t and t' containing /. This choice of notation simplifies the notation 
and avoids unnecessary (but harmless) ambiguities when taking square roots. 

Let C = € C. Note that our conditions imply that C is a (2r)th root 

of unity, and that is a primitive rth root of unity; that is, (C^)* 7 ^ 1 for 
k = 1 ,..., r — 1. For each positive integer i, we define [i] = {C ~ C~*)/(C ~ C~^) 
and, as a special case, [0] = 1. We next define the “bracket factorial” [*]! = 
[i] [i — 1] ... [0]. Note that [r] = 0, and thus [i]! = 0 for all i> r. 

We give every vertex constant weight 


and to each edge e of colour i € I (i.e., for which 0(e) = i) we give the weight 

|e|, = (-l)2*.[2* + l]. 

A triangle / whose three edges have colours i, j,k G I is assigned the weight 

Ifl =('-lF+J+fc + J - k]l ■ [i + k-j]l ■ [j + k-i]\ 

^ ’ [i + j + k + l]\ 

Note that the parity condition and triangle inequalities ensure that the argument 
inside each bracket factorial is a non-negative integer. 

Finally, let t be a tetrahedron with edge colours ioi*ij* 2 ,* 3 )* 4,*5 as indi¬ 
cated in Figure In particular, the four triangles surrounding t have colours 
(*o,ii!* 3 ); (* 0 :* 27 * 4 ), (*i,* 2 ,* 5 ) and ( 13 , 14 , 15 ), and the three pairs of opposite 
edges have colours ( 10 , 45 ), (A, * 4 ) and ( 12 , * 3 )- We define 

T^{t,z) = [z - io - ii - 43]! • [2 - 4 o - 42 - 44]! • [z - 4 i - 42 - 45]! • [z - 43 - 44 - 45] 
K^{t, Z) = [40 -b 4 i -I- 44 -I- 45 - z]! • [40 -b 42 -b 43 -b 45 - z]! ■ [41 -b 42 -b 43 -b 44 - z]! 

for all integers z such that the bracket factorials above all have non-negative 
arguments; equivalently, for all integers z in the range z“ < z < z"*" with 

z“ = max{ 4 o -b 41 -b 43 , 4 o + *2 + * 4 , 4 i + 42 -b 45, 43 -b 44 + *5} ; 
z+ = min{ 4 o -b 4 i -b 44 + 45, 4 o + 42 -b 43 -b 45, 4 i -b 42 -b 43 -b 44}. 


Note that, as before, the parity condition ensures that the argument inside each 
bracket factorial above is an integer. We then declare the weight of tetrahedron 
t to be 


\tu= E 

Z~ + 


T^{t,z) ■ K4>{t,z)' 
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Note that all weights are polynomials on C with rational coefficients, where 
^ '?/»■_ 

Using these weights, we define the weight of the colouring to be 

v^V eeE f^F t^T 

and the Turaev-Viro invariant to be the sum over all admissible colourings 

TV.,,(T) = ^ |T|,. 

9 admissible 

In [19], Turaev and Viro show that TVr,g(T) is indeed an invariant of the 
manifold; that is, if T and T' are generalised triangulations of the same closed 

3- manifold M, then TYr,q[‘^) = TVr,q(T') for all r,q. Although TVr.^q(T) is 
defined on the complex numbers C, it always takes a real value (more precisely, 
it is the square of the modulus of a Witten-Reshetikhin-Turaev invariant) [22] , 

2.2 Treewidth and parameterised complexity 

Throughout this paper we always refer to nodes and arcs of graphs, to clearly 
distinguish these from the vertices and edges of triangulations. 

Robertson and Seymour introduced the concept of the treewidth of a graph 
which now plays a major role in parameterised complexity. Here, we adapt 
this concept to triangulations in a straightforward way. 

Definition. Let T be a generalised triangulation of a 3-manifold, and let T be 
the set of tetrahedra in T. A tree decomposition (X,{Bt}) of T consists of a 
tree X and bags Bt CT for each node t of A, for which: 

• each tetrahedron t G T belongs to some bag Bt\ 

• if a face of some tetrahedron ti G T is identified with a face of some other 
tetrahedron t 2 G T, then there exists a bag B^. with ti,t 2 G B^; 

• for each tetrahedron t G T, the bags containing t correspond to a connected 
subtree of X. 

The width of this tree decomposition is defined as max \ Br \ — 1. The treewidth 
of‘I, denoted tw(T), is the smallest width of any tree decomposition of T. 

The relationship between this definition and the classical graph-theoretical 
notion of treewidth is simple: tw(T) is the treewidth of the dual graph of T, the 

4- valent multigraph whose nodes correspond to tetrahedra of T and whose arcs 
represent pairs of tetrahedron faces that are identified together. 

Figurej^shows the dual graph of a 9-tetrahedra triangulation of a 3-manifold, 
along with a possible tree decomposition. The largest bags have size three, and 
so the width of this tree decomposition is 3 — 1 = 2. 

Definition. A nice tree decomposition of a generalised triangulation T is a tree 
decomposition {X,{Br}) of T whose underlying tree X is rooted, and where: 

• The bag Bp at the root of the tree is empty {Bp is called the root hag)] 
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1,2,4-2,3,4-3,4,5-3,5,6-6,7 


Figure 2: The dual graph and a tree decomposition of a 3-manifold triangulation 


• If a bag Bt has no children, then \Bt\ = 1 (such a Br is called a leaf hag)-, 

• If a bag Bt has two children B,j and B^, then B^. = B,y = B^ (such a B^. 
is called a join bag) ; 

• Every other bag Br has precisely one child Bg-, and either: 

— \Br\ = \Bg\ -|- 1 and Br D Bg (such a Br is called an introduce hag), 
or 

— \Br \ = \Bg\ — 1 and Br C Bg (such a Bg is called a forget bag). 

Given a tree decomposition of a triangulation T of width k and 0{n) bags, 
we can convert this in 0(n) time into a nice tree decomposition of T that also 
has width w and 0 {n) bags [T5] . 

3 Algorithms for computing Turaev-Viro invari¬ 
ants 

All of the algorithms in this paper use exact arithmetic. This is crucial if we 
wish to avoid floating-point numerical instability, since computing TVr,^ may 
involve exponentially many arithmetic operations. 

We briefly describe how this exact arithmetic works. Since all weights in 
the definition of TVj.,^ are rational polynomials in ^ arithmetic 

operations remain within the rational held extension Q(C)- If C is a primitive 
nth root of unity then this held extension is called the nth cyclotomic field. 
This in turn is isomorphic to the polynomial held Q[A]/$„(A), where $„(X) 
is the nth cyclotomic polynomial with degree <f{n) (Euler’s totient function). 
Therefore we can implement exact arithmetic using degree <p(n) polynomials 
over Q. 

If r is odd and q is even, then is a primitive rth root of unity, and 
Q(C) = Q[A]/$i.(A'). Otherwise C is a primitive (2r)th root of unity, and 
Q(C) = Q[A]/$ 2 r(A'). In this paper we give our complexity results in terms of 
arithmetic operations in Q((). 

Let C be an nth root of unity and Q(C) be the nth cyclotomic held. We repre¬ 
sent elements of Q(C) by polynomials of degree at most (p(n), with rational coef- 
hcients, using the isomorphism Q(C) — Q[A]/$„(A). Asymptotically, the Euler 
totient function satishes (p(n) = 0(n). Additions of two polynomials of degree at 
most n are performed in 0(n) operations in Q, and multiplications and divisions 
are performed in 0{M{n)) operations in Q, with M{n) = O(nlognloglogn) [7]. 
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Hence, for fixed r, Turaev-Viro invariants can be computed in 0{N(r) ■ 
r log r log log r) operations in Q using exact arithmetic over cyclotomic fields, 
where N(r) denotes the number of arithmetic operations needed to compute 
TV.,,. 

3.1 The backtracking algorithm for computing TV.,, 

There is a straightforward but slow algorithm to compute TV.,, for arbitrary 
r, q. The core idea is to use a backtracking algorithm to enumerate all admissible 
colourings of edges, and compute and sum their weights. Both major software 
packages that compute Turaev-Viro invariants—the Manifold Recogniser |15j 
and Regina [5]— currently employ optimised variants of this. 

Let T be a S-manifold triangulation, with £ edges Ci,..., . A simple Euler 

characteristic argument gives i = n + v where n is the number of tetrahedra 
and V is the number of vertices in T. Therefore i S 0(n). 

To enumerate colourings, since each edge admits r — 1 possible colours, the 
backtracking algorithm traverses a search tree of 0((r — 1)^) nodes: a node 
at depth i corresponds to a partial colouring of the edges ei,..., e^, and each 
non-leaf node has r — 1 children (one edge per colour). Each leaf of the tree 
represents a (possibly not admissible) colouring of all the edges. At each node 
we maintain the weight of the current partial colouring, and update this weight 
as we traverse the tree. If we reach a leaf whose colouring is admissible, we add 
this weight to our total. 

Lemma 1. If we sort the edges Ci,..., by decreasing degree, the backtracking 
algorithm terminates in 0{(r — 1)^) arithmetic operations in Q(C). 

Proof. The proof is simple. The main complication is to ensure that updating 
the weight of the current partial colouring takes amortised constant time. For 
this we use Chebyshev’s inequality, plus the observation that the average edge 
degree is < 6. 

In more detail, suppose that the edges ei,..., are ordered by decreasing 
degree. Let deg(ei) be the degree of edge e^. Changing the colour of affects 
the colours of the deg(ei) triangles and deg(ei) tetrahedra containing d. Hence 
the update of the current partial colouring weight is performed in 0(deg(ei)) 
arithmetic operations in Q(C). The total number of arithmetic operations per¬ 
formed by the algorithm is consequently 0{'^^{r — I)Meg(ei)). Following an 
Euler characteristic argument, a triangulation of a closed 3-manifold with I 
edges and v vertices has n = i — v tetrahedra and, consequently, the average 
degree of an edge is 6 {£ — v)/£ and thus constant. Considering that the sequence 
((r — I)®)i is increasing and dege^ is decreasing, we conclude using Chebyshev’s 
sum inequality that 

O “ 1)* deg(ei)^ = O = O ((r - 1)^) 

□ 

To obtain a bound in the number of tetrahedra n, we note that a closed and 
connected 3-manifold triangulation with n > 2 tetrahedra must have v < n + 1 
vertices. Combined with n = £ — v above, we have a worst-case running time of 
0((r — l)2"+i) arithmetic operations in Q(C). 
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3.2 A polynomial-time algorithm for r = 3 

Throughout this section, T will denote an n-tetrahedra triangulation of an ori- 
entable 3-manifold M. 

We introduce homology with coefficients in the field Z 2 . A generalised trian¬ 
gulation T, after gluing, contains a set of vertices (dimension 0), edges (dimen¬ 
sion 1), triangles (dimension 2) and tetrahedra (dimension 3) with incidence 
relations. The group of d-chains of T, d G {0,... ,3}, denoted by Cd(T, Z 2 ) is 
the group of formal sums of d-simplices with Z 2 coefficients. The dth boundary 
operator is a linear operator : C£;(T,Z 2 ) —)■ Cd_i(T, Z 2 ) that assigns to a 
d-face the alternate sum of its boundary faces. The kernel of dd^ denoted by 
Zd(T, Z 2 ), is the group of d-cycles and the image of dd, denoted by Brf_i(T, Z 2 ), 
is the group of (d — 1)-boundaries. The fundamental property of homology is 
that dd o dd+i = 0, and we define the dth homology group Hd(T,Z 2 ) of T to 
be the quotient Hd(T, Z 2 ) = Z£;(T, Z 2 )/B(i(T, Z 2 ). It is a finite dimensional 
Z 2 -vector space and we denote its dimension by f3d{^, ^ 2 )- For a more thorough 
introduction into homology theory, see m- 

The value of TV 3 ^g(T), q G {1,2}, is closely related to H 2 (M,Z 2 ), the 2- 
dimensional homology group of M with Z 2 coefficients. H 2 (M, Z 2 ) is a Z 2 -vector 
space whose dimension is the second Betti number /32(M, Z 2 ). Its elements are 
(for our purposes) equivalence classes of 2-cycles, called homology classes, which 
can be represented by 2-dimensional triangulated surfaces S embedded in T. 

The Euler characteristic of a triangulated surface S, denoted by x(5'), is 
x(S') = V — e f, where v, e and / denote the number of vertices, edges and 
triangles of S respectively. We define the Euler characteristic x(c) of a 2-cycle c 
to be the Euler characteristic of the embedded surface it represents. Given T, the 
dimension /32(M, Z 2 ) of H 2 (M, Z 2 ) may be computed in 0(poly(n)) operations. 

The following result is well known M- 

Proposition 2 (Proposition 8.1.7 in [Mj). Let M be a closed orientable 3- 
manifold. Then TV 3 ^ 2 (Af) is equal to the order o/H 2 (M,Z 2 ). Moreover, if 
M contains a 2-cycle with odd Euler characteristic then TV 3 q(M) = 0, and 
otherwise TV 3 q(M) = TY. 

Since H 2 (M,Z 2 ) is a Z 2 -vector space of dimension /32(M,Z 2 ), we have 
TV 3 , 2 (AI) = 2 ^ 2 (m,Z 2 )^ and one can compute TV 3 ^ 2 (AI) in polynomial time. 
The parity of the Euler characteristic of 2-cycles does not change within a ho¬ 
mology class; moreover, if M is orientable, the map H 2 (M,Z 2 ) —)■ Z 2 , taking 
homology classes to the parity of their Euler characteristic, is a homomorphism. 
Consequently, one can check whether TV 3 q(M) = 0 or TV 3 q(M) = TV 3 , 2 (AI) 
by computing the Euler characteristic of a cycle in each of the /32(M, Z 2 ) ho¬ 
mology classes that generate H 2 (M,Z 2 ). Because ,02(^1,Z 2 ) = 0{n), this leads 
to a polynomial time algorithm also. 

3.3 7(/:P-hardness of TV4 1 

The complexity class ffP is a function class that counts accepting paths of a non- 
deterministic Turing machine [20] . Informally, given an NP decision problem C 
asking for the existence of a solution, its ffP analogue ffC is a counting problem 
asking for the number of such solutions. A problem is ffP-hard if every problem 


in polynomially reduces to it. For example, the problem ^3SAT, which 
asks for the number of satisfying assignments of a 3CNF formula, is ^P-hard. 

Naturally, counting problems are “harder” than their decision counterpart, 
and so #P-hard problems are at least as hard as iVP-complete problems— 
specifically, ^P complete problems are as hard as any problem in the polynomial 
hierarchv|18]. Hence proving hardness is a strong complexity statement. 

Kirby and Melvin El prove that computing the Witten-Reshetikhin-Turaev 
invariant is ^P hard for r = 4. This invariant Tr is a more complex 3- 
manifold invariant which is closely linked to the Turaev-Viro invariant TV^p 
by the formula TVr,i(M) = \Tr{M)\^. Although computing TV^^i is “easier” 
than computing t^, the we can adapt the Kirby-Melvin hardness proof to fit our 
purposes. 

To prove their result, Kirby and Melvin reduce the problem of counting the 
zeros of a cubic form to the computation of T 4 . Given a cubic form 

0 ( 3 ^ 1 , . . . , 3^n) ^ ^ G T ^ ^ ^ij A ^ ^ ^ijk XiXjXj^ 

i i,j 

in n variables over Z/2Z and with zeros, they define a triangulation of a 
3-manifold Me with 0(poly(ri)) tetrahedra satisfying Ti{Mc) = 2^c — 2" and 
hence TV^p = ( 2 #c- 2 ") 2 . 

Consequently, counting the zeros of c(a:i,..., Xn) reduces to computing Ti{Mc), 
and so computing TV 4 ^i determines #c up to a ± sign ambiguity (depending 
on whether or not c admits more than half of the input as zeros). 

Establishing the existence of a zero for a cubic form is an AP-complete 
problem, which implies that counting the number of zeros is ^P complete. 
Consequently, computing r 4 is ^P hard. Kirby and Melvin prove this claim 
explicitly by reducing ^3SAT to the problem of counting the zeros of a cubic 
form; moreover, we observe that their construction ensures that this cubic form 
admits more than half of its inputs as zeros. 

We recall the reduction of ff3SAT to the problem of counting the number 
of zeros of a cubic form over Z/2Z found in [TT]. Given a 3CNF formula over 
variables Xi,...,Xn'. 

C = Cl A ... A Cm with Cj = V Xj^ V Xj^ 

and Xi is either Xi or its negation the problem ffZSAT consists in counting 
the number of assignments of ’’true” and ’’false” to the variables a;i,...,a;„ 
satisfying the formula. 

To each form Cj = Xj^ V Xj.^ V Xj^ we assign a cubic equation qj over Z/2Z 
by setting ’’true” = 0 and ’’false” = 1 , and replacing Xi by the variable Xi and 
^Xi by (1 — Xi). For example, a form {-^Xi A —~Xj A Xk) leads to the equation 
(1 — Xi){l — Xj)xk = 0. An assignment satisfies Cj if and only if it cancels 
qj, hence the number of solutions to the system of equations {qi{xi ,..., Xn) = 

0 ,..., qm{xi ,..., Xn) = 0 } is equal to ffc, the number of satisfying assignments 
for C. 

We turn each cubic equation into two quadratic equations by introducing 
a new variable Xij for each monomial XiXjXk of degree 3 and a new quadratic 
equation Xij — XiXj = 0, and by replacing the product XiXjXk by XijXk. We 
obtain a set of m' equations {gi(xi,..., Xn) = 0 ,..., qm'(xi ,..., Xn') = 0 } in 
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n' variables over Z/2Z, with m < m' < 2m and n < n! < 2n. The number of 
solutions of this system remains #c. 

Finally, we define the following cubic form c by introducing m' extra variables 

'^1, ■ • ■ 5 ■ 

m' 

Q — ^ ^ (^1 5 ■ ■ • ) ^n' ) 

i=l 

The number of zeros of Q is equal to 2™ + 2™ “^(2” — ^c) > |2” . 

Because Q is defined on n' + m' variables it admits more than half of its input 
as zeros. Finally, ^3SAT reduces to counting the number of zeros of a cubic 
form which admits at least half of its input as zeros. 

Thus the same reduction process as for T 4 applies for TV 44 , 

Corollary 3. Computing TV 4 P is hard. 

4 A fixed-parameter tractable algorithm 

Here, we present an explicit fixed-parameter algorithm for computing Turaev- 
Viro invariants TVr,q for fixed r. As is common for treewidth-based meth¬ 
ods, the algorithm involves dynamic programming over a tree decomposition 
{X,{Br}). We first describe the data that we compute and store at each bag 
Br, and then give the algorithm itself. 

Our first step is to reorganise the formula for TVr,q(T) to be a product 
over tetrahedra only. This makes it easier to work with “partial colourings” 
corresponding to triangulation edges. 

Definition. Let T be a generalised triangulation of a 3-manifold, and let V, E, 
F and T denote the vertices, edges, triangles and tetrahedra of T respectively. 
For each vertex x € V, each edge x G E and each triangle x S P, we arbitrarily 
choose some tetrahedron A(x) that contains x. 

Now consider the definition of TVr., 5 (T). For each admissible colouring 
9: E ^ I and each tetrahedron t G T, we define the adjusted tetrahedron 
weight \t\'g\ 


\t\e = \t\ex n l^l«x n n \f\o- 

v^V e^E f^F 

A{v)^t A{e)^t A(/)-f 

It follows from equation Q that the full weight of the colouring 9 is just 

m« = ni^r«- 

teT 


Notation. Let A be a rooted tree. For any non-root node r of A, we denote 
the parent node of r by t. For any two nodes cr, r of A, we write u ^ r if cr is 
a descendant node of r. 


Definition. Let T be a generalised triangulation of a 3-manifold, and let V, E, 
F and T denote the vertices, edges, triangles and tetrahedra of T respectively. 
Let (A, {Br}) be a nice tree decomposition of T. For each node t of the rooted 
tree A, we define the following sets: 
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• Tt- C T is the set of all tetrahedra that appear in bags beneath r but not 
in the bag itself. More formally: Tr = (Uo-^t Ba)\Br. 

• Fr C F is the set of all triangles that appear in some tetrahedron t G Tr- 

• Ft C E is the set of all edges that appear in some tetrahedron t GT^.. 

• E* C Er is the set of all edges that appear in some tetrahedron t G Tr 
and also some other tetrahedron t' ^ Tr; we refer to these as the current 
edges at node r. 

We can make the following immediate observations: 

Lemma 4. If t is a leaf of the tree X, then we have Tr = Er = Er = E* = %. 
If T is the root of the tree X, then we have Tr = T, Er = F, Er = E, and 
E* = 0. 

The key idea is, at each node t of the tree, to store explicit colours on the 
“current” edges e G E* and to aggregate over all colours on the “finished” edges 
e G Er\Er. For this we need some further definitions and notation. 

Definition. Again let T be a generalised triangulation of a 3-manifold, and 
let {X, {Br}) be a nice tree decomposition of T. Fix some integer r > 3, and 
consider the set of colours I = {0,1/2,1, 3/2,..., (r — 2)/2} as used in defining 
the Turaev-Viro invariants TVr,g. 

Let T be any node of X. We examine “partial colourings” that only assign 
colours to the edges in Er or E*: 

• Consider any colouring 9: Er ^ I. We call 9 admissible if, for each trian¬ 
gle in Er, the three edges e, f,g bounding the triangle yield an admissible 
triple {9{e),9{f),9{g)). 

• Define T,- to be the set of all colourings tp: Ep ^ I that can be extended 
to any admissible colouring 9: Er ^ T 

• Consider any colouring ip G 'i’r {so ip: Ep ^ I). We define the “partial 
invariant” 

TV,.,(T,r,V')= E n l^l«- 

9 admissible t^T-r 
9—ip on E* 


Essentially, the partial invariant TYr,q (T, t, ip) considers all admissible ways 
9 of extending the colouring ip from the current edges Ep to also include the “fin¬ 
ished” edges in Er, and then sums the partial weights \t\'g for all such extensions 
9 using only the tetrahedra in Tr- 

We can now give our full fixed-parameter tractable algorithm for 

Algorithm 5. Let T be a generalised triangulation of a 3-manifold. We com¬ 
pute TVr,g(T) for given values of r and q as follows. 

Build a nice tree decomposition {X, {Br}) of T. Then work through each 
noder of X from the leaves of X to the root, and compute andTYr,q{T,T,ip) 
for each ip G 'i’r as follows. 

1. If T is a leaf bag, then Ep = Er = 0 , 4',- contains just the trivial colouring 
Ip on 0 , and TYr,q{T,T,ip) = 1. 
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2. If T is some other introduce bag with child node a, then Tr =T^. This 
means that 4',- = and for each tp G we have TVr.^g(T, r, i/') = 

3. If T is a forget bag with child node a, then = T^\J {t} for the unique 
“forgotten” tetrahedron t G B-r\Ba-. Moreover, E* extends E* by including 
the six edges oft (if they were not already present). 

For each colouring ip G 5'<t; enumerate all possible ways of colouring the six 
edges of t that are consistent with ip on any edges of t that already appear 
in E*, and are admissible on the four triangular faces oft. Each such 
colouring on t yields an extension ip': E) ^ I of ip: E* —> I. We include 
Ip' in and record the partial invariant TVj.^g(T, r, ip') = TVr,q(T, <J, ip). 

4 . If T is a join bag with child nodes cri,a 2 , then Tr is the disjoint union 
To-i U Ta-^. Here E) is a subset of E*_^ U E*^ . 

For each pair of colourings ipi G and ip 2 G 'I' 0-2 ; */ 'f’l o,nd ip 2 agree 
on the common edges in E*^ n E*^ then record the pair (ipi,ip 2 )- 

Each such pair yields a “combined colouring” in 5't, which we denote by 
ipi ■ ip 2 ' Ep ^ I; note that different pairs (f/'i, V' 2 ) might yield the same 
colouring ipi • ip 2 since some edges from E*_^ U E*^ might not appear in E). 
Then 'I't consists of all such combined colourings ipi •ip 2 from recorded pairs 
{ipi,ip 2 )- Moreover, for each combined colouring ip G we compute the 
partial invariant TYr,q(fji,T,ip) by aggregating over all duplicates: 

TYr,qi1,T,iP) = TV,.,(T,(Ti,^i)-TV,,,(‘I,a2,V'2). 

recorded 

Once we have processed the entire tree, the root node p of X will have E* = 0 , 
'I'p will contain just the trivial colouring ip on 0, and TYr,q((Z, p, ip) for this 
trivial colouring will be equal to the Turaev-Viro invariant TVr,q(T). 

The analysis of the time complexity of this algorithm is straightforward. 
Each leaf bag or introduce bag can be processed in 0(1) time (of course for the 
introduce bag we must avoid a deep copy of the data at the child node). Each 
forget bag produces < (r — 1)1^* I colourings, each of which takes 0(|E*|) 
time to analyse. 

Naively, each join bag requires us to process |d>cri | ■ \^c! 2 1 < (^ ~ 1 

pairs of colourings (ipi,ip 2 ). However, we can optimise this. Since we are only 
interested in colourings that agree on E*^ n E*^, we can first partition 
and into buckets according to the colours on E*^ n E*^ , and then combine 
pairs from each bucket individually. This reduces our work to processing at 
most (r — 1)1®'^!^^'^ 2 1 pairs overall. Each pair takes 0{\Ep.\) time to process, 
and the preprocessing cost for partitioning is O • log = 

o((r-l)l^-J.|E*j2 logr). 

Suppose that our tree decomposition has width k. At each tree node r, 
every edge in E) must belong to some tetrahedron in the bag B^, and so \Ep\ < 
6 (fc + 1). Likewise, at each join bag described above, every edge in or E*^ 
must belong to some tetrahedron in the bag Hg-. and therefore also the parent 
bag Br, and so \E*_^ U E*^ \ < Q{k + 1). From the discussion above, it follows 
that every bag can be processed in time O ((r — l)®(^+i) . k^ logr), and so: 
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Theorem 6. Given a generalised triangulation T of a 3-manifold with n tetra- 
hedra, and a nice tree decomposition of T with width k and 0{n) bags, Algo¬ 
rithm computes TVr,q(T) in O {n ■ (r — • k^ logr) arithmetic opera¬ 

tions in Q(C)- 

Theorem shows that, for fixed r, if we can keep the treewidth small then 
computing TV^^g becomes linear time, even for large inputs. This of course is 
the main benefit of fixed-parameter tractability. In our setting, however, we 
have an added advantage: TVr,q is a topological invariant, and does not depend 
on our particular choice of triangulation. 

Therefore, if we are faced with a large treewidth triangulation, we can retri¬ 
angulate the manifold (for instance, using bistellar flips and related local moves), 
in an attempt to make the treewidth smaller. This is extremely effective in prac¬ 
tice, as seen in Section 

Even if the treewidth is large, every tree node r satishes \E*\ < i, where I 
is the number of edges in the triangulation. Therefore the time complexity of 
Algorithmj^reduces to O (n • (r — 1)^ • log r), which is only a little slower than 
the backtracking algorithm (Lemma [^. This is in sharp contrast to many FPT 
algorithms from the literature, which—although fast for small parameters— 
suffer from extremely poor performance when the parameter becomes large. 


5 Implementation and experimentation 

Here we implement Algorithm]^ (the fixed-parameter tractable algorithm), and 
subject both it and the backtracking algorithm to exhaustive experimentation. 

The FPT algorithm is implemented in the open-source software package 
Regina [5]: the source code is available from Regina’s public git repository, and 
will be included in the next release. For consistency we compare it to Regina’s 
long-standing implementation of the backtracking algorithmj^ 

In our implementation, we do not compute treewidths precisely (an NP- 
complete problem)—instead, we implement the quadratic-time GreedyFillln 
heuristic [2], which is reported to produce small widths in practice [21] . This 
way, costs of building tree decompositions are insignificant (but included in the 
running times). For both algorithms, we use relatively naive implementations of 
arithmetic in cyclotomic fields—these are asymptotically slower than described 
in Section]^ but have very small constants. 

We use two data sets for our experiments, both taken from large “census 
databases” of 3-manifolds to ensure that the experiments are comprehensive 
and not cherry-picked. 

The first census contains all 13400 closed prime orientable manifolds that 
can be formed from n < 11 tetrahedra H H]. This simulates “real-world” 
computation—the Turaev-Viro invariants were used to build this census. Since 
the census includes all minimal triangulations of these manifolds, we choose the 
representative whose heuristic tree decomposition has smallest width (since we 
are allowed to retriangulate). 

The second data set contains the first 500 (much larger) triangulations from 

^The Manifold Recogniser ca also implements a backtracking algorithm, but it is not 
open-source and so comparisons are more difficult. 
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Figure 3: Running times for the 13400 closed census manifolds 


the Hodgson-Weeks census of closed hyperbolic manifolds [8]. This shows per¬ 
formance on larger triangulations, with n ranging from 9 to 20. 

Figures and compare the performance of both algorithms for each data 
set. All running times are for TVyp (the largest r for which the experiments 
were feasible), and are measured on a single 3 GHz Intel Core i7 CPU. Both plots 
use a log-log scale with one data point per input triangulation. The results are 
striking: the FPT algorithm runs faster in over 99% of cases, including most of 
the cases with largest treewidth. In the worst example the FPT algorithm runs 
3.7 X slower than the backtracking, but both data sets have examples that run 
> 440 X faster. It is also pleasing to see a clear impact of the treewidth on the 
performance of the FPT algorithm, as one would expect. 


6 An alternate geometric interpretation 


In this section, we give a geometric interpretation of admissible colourings on a 
triangulation of a 3-manifold T in terms of normal arcs, i.e., straight lines in the 
interior of a triangle which are pairwise disjoint and meet the edges of a triangle, 
but not its vertices (see Figure]^. More precisely, we have the following 

Theorem 7. Given a 3-manifold triangulation T, and r > 3, an admissible 
eolouring of the edges of T with r — I colours corresponds to a system of normal 
arcs in the 2-skeleton with < r —2 arcs per triangle forming a collection of cycles 
on the boundary of each tetrahedron of T. 


Proof. Following the definition of an admissible colouring from Section 2.1 the 
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Figure 4: Running times for the first 500 Hodgson-Weeks census manifolds 


colours of the edges ei, 62 , 63 of a triangle / of 1 must satisfy the parity 
condition, the triangle inequalities, and the upper bound constraint. 

For a colouring 9(e) of an edge e, we define ^(e) = 29(e) which is an in¬ 
teger; we also use the term ’’colouring” for (j). We interpret the colourings 
0(ei), 0 ( 62 ), (/>{e 3 ) G {0,1,..., r — 2} as the number of intersections of normal 
arcs with the respective edges of the triangulation (see Figure [^. Without loss 
of generality, let (j)(si) > ^( 62 ) > We construct a system of normal arcs 

by first drawing 1 ^( 62 ) arcs between edge Ci and 63 and (l)(^i)~<P(£^ 2 ) arcs between 
edge ei and 63 . This is always possible since (()(ei) < <('( 62 ) +<^( 63 ) by the trian¬ 
gle inequality. Furthermore, the parity condition ensures that an even number 
of unmatched intersections remains which, by construction, all have to be on 
edge 63 . If this number is zero we are done. Otherwise we start replacing normal 
arcs between ci and 62 by pairs of normal arcs, one between ei and 63 and one 
between 62 and 63 (see Figure]^. In each step, the number of unmatched inter¬ 
section points decreases by two. By the assumption 4>(e2) > <(>( 63 ), this yields a 
system of normal arcs in / which leaves no intersection on the boundary edges 
unmatched. This system of normal arcs is unique for each admissible triple of 
colours. By the upper bound constraint, we get at most r — 2 normal arcs on /. 

Looking at the boundary of a tetrahedron t of T these normal arcs form a 
collection of closed cycles. To see this, note that each intersection point of a 
normal arc in a triangle with an edge is part of exactly one normal arc in that 
triangle and that there are exactly two triangles sharing a given edge. □ 

Now, let T be a closed n-tetrahedron 3-manifold triangulation, t a tetrahe¬ 
dron of T, /i and /2 two triangles of t with common edge e of colour 4>(e), and 
Oi and bi the respective non-negative numbers of the two normal arc types in fi 
meeting e, i G {1,2}. Since the system of normal arcs on t forms a collection of 
cycles on the boundary of t, we must have oi-I -61 = 02 - 1-62 < '^—2, giving rise to 
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(j>{ei) = 5 ; 


0(63) = 3 


Figure 5: Constructing a system of normal arcs from edge colourings. 


a total of 6n linear equations and 12n linear inequalities on 6n variables which 
all admissible colourings on T must satisfy. Thus, finding admissible colourings 
on T translates to the enumeration of integer lattice points within the polytope 
defined by the above equalities and inequalities. 

Now, if we drop the upper bound constraint above, we get a cone. Computing 
the Hilbert basis of integer lattice points of this cone yields a finite description 
of all admissible colourings for any r > 3 and, thus, the essential information to 
compute TVr,q(T) for arbitrary r. Transforming this approach into a practical 
algorithm is work in progress. 
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